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Abstract 

We study the stationary state of a rough granular sphere immersed in a thermal 
bath composed of point particles. When the center of mass of the sphere is fixed 
the stationary angular velocity distribution is shown to be Gaussian with an effective 
temperature lower than that of the bath. For a freely moving rough sphere coupled 
to the thermostat via inelastic collisions we find a condition under which the joint 
distribution of the translational and rotational velocities is a product of Gaussian 
distributions with the same effective temperature. In this rather unexpected case we 
derive a formula for the stationary energy flow from the thermostat to the sphere in 
accordance with Fourier law. 
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1 Introduction 



The study of the dynamics of a tracer particle is a classical problem of nonequilibrium 
statistical mechanics. The object is to determine the evolution of the state of a single 
particle resulting from interaction with the surrounding medium. Such a study can yield 
precious information on the effects of many-body dynamics. Moreover, the relative sim- 
plicity of the problem creates chances for precise theoretical predictions. It has recently 
motivated a number of works aiming at understanding the dynamics of fluidized granular 
media |lj. In particular, the evolution of a granular sphere immersed in a granular medium 
homogeneously cooling down has been discussed j2] as well as the Brownian motion in a 
granular fluid [3J. The case of an impurity put in a vibrated low-density granular system 
has also been studied 0]. 

An interesting qualitative question related to the effects of inelastic collisions taking 
place in granular fluids is that of the resulting structure of the distribution function when 
different kinds of degrees of freedom are present. A gas of rough spheres with both trans- 
lational and rotational degrees of freedom has been recently examined from this point of 
view [5]. The main prediction based on numerical studies and approximate analytic argu- 
ments is that dissipative collisions induce statistical dependence between orientations of 
the angular and translational velocities. 

In the present paper we address an analogous question in an even simpler situation of a 
single tracer granular rough sphere suffering inelastic collisions with point masses forming a 
low-density thermal bath. Our object is to find out what kind of stationary state can result 
from a dissipative coupling to a thermostat. In the case of elastic collisions the particle 
would eventually attain equilibrium at the temperature of the bath. But the nature of 
the asymptotic stationary state in which there is a constant dissipative heat flow from the 
thermostat to the tracer particle remains a largely open question of fundamental interest. 

Exact results have been derived for a smooth hard sphere where inelastic collisions 
could influence only translational motion. It turned out that at the level of the Boltzmann 
kinetic theory the stationary velocity distribution had the form of a Maxwell distribution 
with an effective temperature lower than that of the thermostat j6j. In one dimension 
one could even rigorously solve the initial value problem deriving in particular the exact 
dependence of the diffusion coefficient on dissipation [7j. In the case of purely translational 
degrees of freedom the appearance of a Gaussian distribution has been shown to follow from 
the equivalence between the Boltzmann equation for a granular tracer particle suffering 
inelastic collisions and the Boltzmann equation for an elastic tracer particle with a suitably 
modified mass [2]. 

This remarkable property also occurs if the test particle is rough but has a fixed mass 
center and thus only rotational degrees of freedom [8] (an original derivation of this fact 
is presented in Section 3). The stationary angular velocity is then again Gaussian with an 
effective temperature lower than that of the thermostat. 

In the present paper we extend the study of stationary states at the level of Boltzmann's 
kinetic theory to the case of a tracer rough sphere whose both translational and rotational 
motions are influenced by inelastic collisions (the distribution of kinetic energy in a granular 
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gas composed of rough spheres has been discussed in p|, [10J ). In section 2 we describe 
the model. Section 3 contains the description of our method first illustrated on simple 
situations where only one type of degrees of freedom is present. We then turn to the general 
case and show that when the restitution coefficients for the normal and tangential relative 
velocities obey a specific relation (I3T1) the joint velocity distribution becomes a product of 
two Maxwell distributions for the angular and translational velocities corresponding to the 
same effective temperature. It turns out that the derived relation may be fulfilled only if 
the restitution coefficient relevant for rotational motion is larger than that linked to the 
motion of the mass center. 

The occurence of a stationary factorized Gaussian distribution for the two types of 
degrees of freedom inelastically excited by collisions is quite remarkable, and, in view of 
the results obtained for a gas of granular rough spheres [5] rather unexpected. The main 
result of Section 3 is obtained by using an appropriate change of integration variables in 
the gain term of the Boltzmann equation ( the method generalizes that used in Ref. [2] ). 

The heat flux that maintains the test particle in a stationary state is calculated in 
Section 4. It obeys the analogue of Fourier law with a thermal conductivity proportional to 
the temperature jump between the sphere and the thermostat, as in generic hydrodynamic 
theories. In Section 5 we briefly comment on the possible structure of the joint velocity 
distribution when the values of the restitution coefficients are not related by the equation 
derived in Section 3. We expect that the typical case would involve statistical dependence 
between the angular and translational velocities of the sphere. 



For the sake of simplicity the thermal bath particles are supposed to be point masses m 
performing purely translational motion. Their distribution in the one-particle phase space 
is the product of a uniform spatial density p and a Maxwell velocity distribution 0t( v ; m ) 
corresponding to temperature T 



D is the dimension of the space (D = 2 or 3), and ks is Boltzmann's constant. 

The rough sphere is supposed to have mass M, radius R, moment of inertia J, and to 
move with translational velocity V, and angular velocity O. Its total kinetic energy equals 
thus 



where q = I /MR 2 is a number reflecting the mass density distribution inside the sphere 
(disk). 



2 The model 




(i) 




(2) 
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2.1 Collisional laws 



Consider a binary collision between the rough sphere and a point particle of the thermostat. 
The instantaneous collisional transformation of velocities 



(V,n,v)-(V*,n*,v*) (3) 

is conveniently described with the help of the unit vector n along the line segment from 
the center of the sphere to the point of impact. 

The linear velocity of the point at the surface of the sphere hit by the thermostat particle 
is (V + Rfl x n) . The relative velocity at which the particle approaches the impact point 
is thus 

g = v- V-ffixE (4) 
In what follows we will use the notations 



A n = (A • n) n, and A t = A — A n = n x (A x n) (5) 

for the normal and tangential components of any vector A. 

In the simplest model of inelastic collisions one defines the instantaneous transformation 
of relative velocity (f4]) by (see e.g.pQ) 

g^ = -agn (6) 

St = -P&t (7) 

where < a < 1, and — 1 < /3 < 1 are the translational and rotational restitution 
coefficients, respectively. The cases of a = 1 with /3 — ±1 correspond to elastic encounters. 

In order to calculate postcollisional velocities (V*,0*,v*) one has to use relations 
(E]),(I7j) together with the conservation laws for both the momentum and the angular mo- 
mentum. One then finds (see e.g. [10J) 

RSl ' = m+ i + 1 ,a + M) iixs ' 

V = v + i±^ g „ + -^l±|) gl (8) 
1 + /i 1 + g(l + n) 



l + l + g(l + /i)° 

where fi denotes the mass ratio fi = M/m. 

According to the defining equations (EJ) and (CO) the formulae corresponding to the 
inverse collision 

(v**, ir*, v**) -> (v,n,v) (9) 

are obtained by changing a into 1/a, and f3 into 1/(3. 
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2.2 The Boltzmann equation 

We denote by /(V, SI, t) the probability density for finding the test particle at time t with 
translational velocity V and angular velocity SI. 

Within Boltzmann's theory whose predictions are the object of the present paper, the 
tracer particle never suffers recollisions, and thus always encounters unperturbed thermal- 
ized particles. The Boltzmann kinetic equation is thus linear in this case and in homoge- 
neous situations it has the form 

dt = pR 



Jdnjdve([v- V] -n)| [v- V] -n| (10) 



a ,, 

In the gain term there appears the factor l/a 2 \f3\ D ~ 1 which guarantees the conservation of 
the normalization of the velocity distribution. Q(x) denotes here the unit step function. 
We restrict further discussion to the stationary solution of equation (|10|) . 

It will be useful to display the detailed structure of the gain term. To this end we 
express the precollisional velocities (V**, ft**, v**) in terms of the scaled relative velocities 

Wn = = i (V - V)„ (11) 

w t = ig t = ^ (v - V - RSI x n) t (12) 
Using equations (jHJ) with a and (3 replaced by or 1 and we find 

V- = V + i±gw. + f +f w. (13) 
1 + fi l + g(l + /i) 

v ** = v + flnxn-w + — w n + \ — ; V w * 

1 + /i 1 + g(l + /i) 

Under the change of the integration variables (v n , v t ) — > (w n , w t ) the gain term takes the 
form 

/poo P 
dh / dw n w n / dw t (14) 



1 + ^ l + g(l + /i) r ' + g(l + //)] 

X0 T V + RSI x n - w + — — w n + \ Hi \ f J w t ; m 
V 1 + 1 + g(l + n) 

where w n = \w n \, and the relation dv n dvt = a|/5| £) ~ 1 (iw ri (iw i has been taken into account. 
The use of the scaled relative velocities (ITT]) . ( 1T21) as integration variables will play a crucial 
role in exploring the possibility of finding an elastic problem equivalent to the dissipative 
one. 
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2.3 Energy dissipated in a binary collision 

The collisional transformation laws (jHJ) conserve the center of mass velocity 

MV + mv /jY + v 



V 



° M + m 11+ 1 1 

the tangential vector 

C t = v t - V t + g(l + x n, (15) 

and also the normal component of the angular velocity Q n . Moreover, the transformation 
laws for the relative velocity g reduce to simple rescalings (jEJ), (ED). 

The evaluation of the collisional change in the total kinetic energy of the colliding pair 

AE tot = E tot (v\ n*, v*) - E tot (v, n, v) 

where 



E tot (V, Q, v) = E(V, n) + -mv 2 = -m j/iV 2 + -g/xi? 2 ft 2 + v 2 j , (16) 

can be thus conveniently performed by expressing E tot in terms of the invariant vectors 
Vg, f^n, Ct, and the relative velocity g. Owing to the relations 

v 2 + /iV 2 = (1 + /x) V 2 , + ( V _ V) 2 (17) 



and 



we find 



(v - V) 2 + q(l + /U ) J R 2 O t 2 = 1 + + [q(l + /i)g 2 + C 2 ] (18) 



E tot = { + [g 2 + ,(1 + ,)^ 2 ] + YT ^ ) W + ,)g 2 + C 2 ] } 

(19) 

The dissipation of the total energy E tot under binary collisions is due to the reduction of 
the length of vectors g n and g t governed by the restitution coefficients a and (3. In view 
of the invariance of vectors Vg, f2 n , and Ct one finds 

m V f ft 2\ „2 , + /, ,o2x 2 



A£ « - " t (Tm C 1 " + rfioTT) (1 " " ^ > (20) 

The total pair energy is conserved if a = 1 and /3 = ±1. 



6 



3 Solvable situations 



3.1 Specific elastic cases 

In the case of an elastic rough sphere a = (3 = 1, the relative velocity g is reversed as 
the result of the impact, and the conservation of the kinetic energy AE tot = implies 
the equilibrium solution of the Boltzmann equation involving the product of two Gaussian 
distributions 

/ e «(V,fi) = «MV;M)0™*(fi;I) (21) 
satisfying the equipartition law. Here 

/ i \ D/2 ( m 2 \ 

exp (-^) (22) 

with / defined in Eq.([2]). The Maxwell distribution 0r(V; M) has been defined in Eq.([T]). 

The case of a — 1, /3 — — 1 corresponds to an elastic smooth sphere: the rotational 
degrees of freedom are not influenced by collisions. Consequently, any distribution of the 
form 

T (V;M) X (ft) (23) 
represents a stationary state of a smooth sphere. 



3.2 Connection with an elastic problem for a rough sphere 

Let us consider now the possibility for reducing the Boltzmann equation for a granular 
sphere (inelastic collisions) to an equivalent case of a rough sphere elastically coupled to 
the thermostat. The main observation is that the gain term ( fl4l) would take the elastic 
form for a rough sphere (a = 1 and (3 = 1) with a mass Jim 

PR " 1 [ dh [ d Wn W n [ dwj ( V + -^W„ + - ; ^ ; zr W t , fl ^ X " ' 



1 + l + g(l + /2) " R[l + q(l + JZ)} 

x0 T ( V + Ml x n - w + -^-cw w + 2 ^ | ^ w t ; m) (24) 
V l + H l + q(l + ii) J 

provided the equations 

I±^ = -^= and 1 + ^ 2 f2 5) 

l + 1 + + 1 + g(l + /x) 

could be simultaneously satisfied. This simple fact is the basis for the further discussion. 
The corresponding formula for the loss term 

/r+oo /• 
rfn J dg n g n J dg t /(V, ft)0 T (V + RQ x n + g; m) (26) 
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is irrelevant from the point of view of the correspondence with an elastic case as it does 
not involve the collision law. 

Before discussing the general case let us consider two simple situations where only one 
type of degrees of freedom is present. 

When /3 — — 1, the sphere is smooth and only the translational velocity V is influenced 
by collisions. We are thus interested in the stationary velocity distribution / s t(V). For a 
comparison with a smooth sphere undergoing elastic collisions with the point particles of 
the thermostat, we make the extra change of variable w t — > —w t in the gain term (|14p 
which becomes for a smooth sphere {(3 = —1) 



G(V, ft) = x(fl)pR D ~ l I dn I dw n w n I dw t (27) 



oo 







/ ( v + y^ w «) <fa ( v + Rft x A + w - ( 2 - ) w «; m 



(The complete variable change is similar to that used in Ref. [2].) The gain term (12T 
would take the elastic form for a smooth sphere (a = 1 and (3 = — 1) with a mass Jim 



OC 



x(fl)pR J dhj dw n w n J dw t (2S 

f(y+ T^ w «) fa ( V + m x fi + w - T~^ w " ; m 

provided that (1 + a)/(l + /x) = 2/ (1 + ju). Eventually, from the two conditions ( |25l) only 
the first remains relevant yielding the mass ratio 

p = P+ l —^{l + p) (29) 
1 + a 

The stationary distribution of the mass center velocity has thus the form of the Maxwell 
distribution at temperature T for a particle with mass M = Jim, or, equivalently, it is a 
Maxwell distribution for the sphere of mass M with an effective temperature T* satisfying 
the relation M/T = M/T*. From §M) we find 

T,= (1 + a) ? T < T (30) 
1 — a + 2/i 

Hence / s t(V) = 0r*(V;M) which is the result derived in [6]. 

When the rough sphere has a clamped center of mass, V = 0, only rotational degrees 
of freedom are involved. The corresponding collision laws can be deduced from equations 
(jHJ) by taking the limit fi — > oo and g — > 0, while keeping the product g/i = I/mR 2 fixed. 
Indeed, in this limit the moment of inertia is kept unchanged while the mass relevant to 
translational motion tends to infinity which immobilizes the center of mass of the sphere. 
We are thus interested in the stationary velocity distribution f st (fl). The gain term ([TJ 
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takes the form 



G(fi) = pR u ~ L J dhj dw n w n J dw t (31) 
/ (ft — — Wf x n j 0^ ( i?0 x n — w H — w t ; m 



R[l + qiA J \ 1 + qfJ- 

The gain term (13T1) would take the elastic form for a clamped rough sphere (V = 0, a = 1 
and /? = 1) with a mass /Im 

pR ' 1 dh dw n w n / c?w t 



x/ H- WjXnUr fflxn-w+ ~ w f ;m (32) 

provided that 

1 + 2 

- — = = 33 

1 + qfi 1 + q\x 

Eventually in the above mentioned limit only the second of the two relations (125jl remains 
relevant : it takes the limit form ( 1331) and yields the mass ratio 

_ 1 - (3 + 2g/x 

^TTT^T^ (34) 

The stationary solution of the Boltzmann equation has thus again the form of a Maxwell 
distribution for the angular velocity f2 



( q Mg*\ D > 2 



with the effective temperature given by 



(35) 



2k B T^ j 



o J1+^ T < T (36) 
** 1 - /? + 2q/M v 7 

The results ( 135]) . ( 1361) have been previously derived in an unpublished work [8] in a different 
way. We notice that is independent of a, whereas the collisional laws are not. 

In the general case of a granular rough sphere (—1 < (3 < +1, < a < 1) with 
translational and rotational degrees of freedom both conditions (1251) are to be satisfied 
simultaneously implying the relation between the restitution coefficients of the form 

P = a + ] —^ 7 (37) 

The relation can occur only for positive values of the restitution coefficient (3 that are larger 
than a. It is thus possible to rewrite the gain term (1141) in the form corresponding to elastic 
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collisions provided equation (137)) is satisfied. This fact is quite remarkable, and, as we have 
mentioned in the introduction, rather unexpected. It would be interesting to understand 
why the particular tuning ()37j) of the parameters governing the dissipation makes appear 
the Gaussian distribution. 

We notice that the constraint ( 1371) can be rewritten as T* = T** with T* given in (1301) 
and T** defined by 



T 



(3* 



(l + g)(l-/3) + 2 9 /i 

The notations have been chosen in order to take into account the property 

T^l = lim g _^o,<?At fixed (see (l36l) ). When TU is expressed in terms of T^, the condition 

XI = TL takes the form 



T. t 



(3 Tl 



-l 



. 1 + MIWJ T " (39) 

where each temperature T* or depends only on one restitution coefficient, a or (3 
respectively. In other words, the constraint ( )37j) expresses the relation ( 139)) between the 
effective temperatures found in the study of a sphere for which collisions could influence 
either only translational or only rotational motion. Our analysis shows that the gain term 
in the Boltmann equation can be rewritten in the form corresponding to elastic collisions 
only when the dissipation related to pure translational motion and the dissipation related to 
pure rotations would lead to effective temperatures linked by equation (139)) . The stationary 
state of the granular rough sphere under the condition (l37j) complies thus with the energy 
equipartition and reads 



M \ D/2 f qMR 2 




k B T, 



eff 



^MV 2 + ^qMR 2 n 2 



(40) 



/ M 
with T eff = T* = T**. 

Clearly, the equivalence with the elastic problem on curve ( )37j) also holds for the Boltz- 
mann equation in the presence of an arbitrary inhomogeneous nonequilibrium state of 
the medium surrounding the test particle, since the structure of the collision term in this 
general case remains unchanged. It seems worth mentionning here that also the Enskog 
equation could be reduced to its elastic form on curve (1371) . since it involves exactly the 
same collisional transformation of velocities as the Boltzmann equation (see e.g. [2]). 
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4 Analogue of Fourier law in the solvable stationary 
state 

In the stationary state determined by the relation ( 1371) the mean kinetic energy of the 
test particle results from the energy equipartition between the rotational and translational 
degrees of freedom 



dv / dnf st (v,n) 



-MV 2 + -qMR 2 fl 2 
2 2 



Dk B T eff (41) 



(D is the dimension of the position space.) This stationary state is sustained by an energy 
flux from the heat bath to the test particle that compensates the energy dissipation under 
binary collisions. 

The collisional rate of the energy dissipation equals 

Jf" = J dV J dn J dv J dnQ ([V - v] ■ n) [V - v] ■ n AE tot f st (V, f2)0 r (v; m) (42) 

Inserting here AE to t calculated in equation (l20l) one finds 



Jfs = I 0^ (2(, + 9 eff ) (D-l)[ g „ + (l + g )fl.„] _ 1 

E 2y /i\l + /i l ; l + q + q/i 1 H) ) 

I 2irk B T . . 

xk B Td — 43 

V m 

where 9 eff = T eff /T and T eff is defined in (HP]). Since a and (3 obey the constraint fl37j) . 
the relations flHUl) and flHHj) allow one to express a and /3 in terms of 6 eff . In particular 

^'iS&v- '■>>»■'> < 44 > 

and 

- 4 [^ + (l + g)M 2(1 "^ /) ^ (45) 
Eventually one finds that the heat flux obeys the analogue of Fourier law 

Jf = k(T) (T e// - T) (46) 

in which the discontinuous jump of the temperature plays the role of the temperature 
gradient present in the theory of continuous media. The thermal conductivity for D = 2 
or 3 reads 



where 9 eff is a dimensionless function of a and // given by ( |30l) or equivalently a function 
of /3 and g/x as given by 
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5 Open questions 



Outside the curve ( 1371) in the (a,/?) plane the possibility of deriving the analytic form of 
the stationary velocity distribution remains an interesting open problem. The simplest 
generalization of solution (HUj) would be the product of two Gaussian distributions with 
different effective temperatures for translational and rotational motions. However, on the 
basis of detailed calculations we conjecture that such a distribution is not possible, and 
that the Gaussian form can appear only on the derived curve ([37j). The complete proof of 
the conjecture is still to be constructed. 

Since the translational velocity V is a vector whereas the angular velocity Q is a 
pseudovector, the stationary state can depend in general on three scalar variables |V|, |f2| 
and |V • 0|, where V and fl are the unit vectors in the direction of V and ft respectively. 
We conjecture that when the constraint (1371) is not respected the stationary distribution 
depends on the angle between V and f2 via the variable |V - fi| which introduces statistical 
dependence between the two velocities. 
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